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Habitat with its components is the most important need of wildlife species. So that, 
habitat management is a critical procedure in wildlife conservation. In the present study, 
desirable habitats of Kavir National Park for Wild goat are determined in two 
distinguished warm and cold seasons, applying Random Forests, Maximum Entropy, and 
Classification Tree algorithmsa and for this matter “Salford Predictive Modeler” and 
MaxEnt software were used. Consequently, an ensemble habitat suitability map was 
provided based on weight-averaging method. The results illustrated that elevation, water 
resources, vegetation type, and slope are the most important factors in predicting Wild 
goat’s suitable habitats with a slightly different impact on each season. In general, among 
the existing plant species, Wild goat prefers Artemisia sieberi, Stipa arabica, Zyghophylom 
eurpterum, and Sieditzia rozmarinus. Also, the altitudes higher than 1050 meters above sea 
level, slopes over 20 to 35 percent, and distances less than 5,000 meters to water 
resources ate desirable. Meanwhile, several suitable patches for Wild goats were found 
which are greatly in need of conservation. Considering the importance of water resources 
in this region, in order to improve the quality of the habitat, further study on water 
resources status is suggested. 
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INTRODUCTION 

Wild goat (Capra aegagrus), taxonomically is included in Order Artiodactyla and Family Bovidae. A 
study by Nadert e¢ a/. (2008) demonstrated that Capra aegagrus is one of the first artiodactyl species 
being domesticated in more than 10,000 years ago in the Fertile Crescent and is the ancestor of 
domestic goat. This species is of high economic and ecological importance, meanwhile, is listed as 
vulnerable in the IUCN Red List (Weinberg, 2008). This species is considered an important prey for 
leopard and due to the status of leopard, being endangered, assessing the Wild goat's habitat is 
crucial for management and conservation processes in different views. 

Today researchers use different species distribution models to conclude the presence and 
location of different species. The models can be the basis of species conservation and decision- 
making for regional development (Fook ef a/., 2009). 

One of the main reasons of wildlife being on the wane, is habitat destruction, in a way that 
many species are forced to live in isolated and semi-isolated populations in fragmented patches 
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(Hirzel & Guisan, 2002). Hence, it is advisable to determine suitable habitats and to act accordingly. 
Nonetheless, in some cases, especially when encountering large-scale regions, there are some 
difficulties such as limited time and narrow budgets that make such studies difficult, if not 
impossible, to take place. Therefore, scientists developed some methods to facilitate such studies, 
and currently, these methods are widely being used in wildlife management (Peterson & Vargas, 
1993). These types of models examine the quantitative and qualitative relationship between the 
ptesence of species and environmental variables because the quality of a habitat results from 
interactions between biological and non-biological variables (Kocev ef a/., 2010). 

In various regions, the crucial factors do not affect habitat suitability for Wild goat, similarly. 
Several studies on Wild goat’s suitable habitats in Iran were carried out using logistic regression, 
Ecological Niche Factor Analysis (ENFA), and Maximum Entropy (MaxEnt) and showed that 
topographic variables, water, and human resources, and vegetation cover were the most important 
factors for habitat selection by a Wild goat (Hosseini ef a/., 2016; Naderi ef a/., 2014; Sarhangzadeh e¢ 
al., 2013; Shams Esfandabad ef a/., 2010). In a study in alpine habitats, Colorado, habitat models of 
mountain goat (Oreamnos americanus) were developed applying the logistic regression method. The 
results demonstrated that mountain goats prefer areas on moderate slopes, near escape terrain, at 
middle-elevations, and on southerly exposures (Gross ef a/., 2002). In another study, in southwestern 
coastal British Columbia, the results indicated that both sexes, especially males, mostly prefer mature 
and old forests (Taylor e¢ a/., 2006). 

The purpose of this research is to evaluate the habitat suitability of Capra aegagrus in Kavir 
National Park, based on the presence data of this species. We categorized a year into two warm (the 
first six months of the year) and cold (the second months of the year) seasons. We attempted to 
ascertain the most essential parameters and to highlight the suitable areas for the species in Kavir 
National Park for the first time, by using the ensemble modeling technique. It is noteworthy that 
neither “Salford Predictive Modeler” (SPM) software nor an ensemble technique has ever been used 
to determine this species habitat in Iran. In the end, some suggestions are presented for 
conservation purposes. 


MATERIAL AND METHODS 

Study area 

This study has been done in Kavir National Park, located in central plateau of Iran (Fig. 1). Kavir, 
with an area of about 4200 square kilometers, is the second largest National Park in Iran, after 
Urmia Lake (Atarodi, 2010). This area is located in latitude 51 degrees 25 minutes to 53 degrees and 
5 minutes E, longitude 34 degrees and 20 minutes to 35 degrees and 10 minutes N. 


Habitat suitability modelling 

In this study, three methods based on machine learning have been used. Machine learning, especially 
predictive models, is increasingly being used in the structure of ecological models (Dzeroski, 2001). 
The three methods used in this project are Classification Tree Algorithms (CTA) (Breiman ¢f a/, 
1984), Random Forest (RF) (Breiman, 2001), and MaxEnt (Phillips e¢ a/, 2006; Elith e¢ a/. 2010). 


Sampling and field studies 

Collecting data started from August 2017 and ended in July 2018. Indices such as horn, fur, teeth, 
feces, carcasses, etc. as well as direct observations of the species were recorded by Global 
Positioning System (GPS), as points of presence. Sampling was carried out in different 
environmental strata of the region (Fig. 2). 
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FIGURE 2. Presence records of Wild Goat in Kavir National Park. 


Bias in datasets can have a significant influence on results. If modeling is performed with bias, it will 
be better to be careful when interpreting the results (Leitao ef a/, 2011). There are two common 
ways to reduce the impact of bias: 1) refining some of the points of presence and 2) defining the 
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background points (pseudo-absence points). The results indicated that the former is more 
appropriate if the number of dataset is big enough to remove some data (Kramer-Schadt, 2013). 
Hence, in present study, the presence points were refined through a spatial filtration by keeping one 
occurrence point within a radius of 1 kilometer. This process was done via "Spatially Rarify 
Occurrence Data" command in SDMtoolbox, ArcMap 10.5 software (Brown, 2014). 


Sample size 

In a study, it was shown that in machine learning models, with only 10 sample points, the accuracy 
of the model, was at 90% of its highest accuracy, and with 50 samples it reaches near to the highest 
accuracy (Stockwell & Peterson, 2002). Figure 2 illustrates the refined presence points in both warm 
and cold seasons. The numbers of presence records in the region are 52 in the warm season (the six 
warm months) and 46 in the cold season (the six cold months). 


Habitat variables 

According to behavioral and ecological characteristics of Capra aegagrus, previous literature 
and the characteristics of the study area, the potential variables in affecting Wild goats distribution 
were chosen as follows: topographical variables including altitude, slope, and aspect, vegetation 
cover variables (density and type), distance to water resources, anthropological variables such as 
distance to roads and DoE guard stations, and climate variables containing monthly minimum and 
maximum average temperature, each for cold and warm seasons respectively, and monthly average 
rainfall. 

For this purpose, some data were collated from Kavir National Park Bureau and 
Meteorological Organization. There are no meteorological and weather stations in Kavir National 
Park, so data from six stations outside the region were interpolated (Childs, 2004) to obtain 
continuous values of climatic variables inside the park. Interpolation was carried out using the ‘IDW’ 
method in ArcMap version 10.5 and climatic layers were produced. According to Khosravi et a/. 
(2014), the optimum interpolation methods for obtaining average temperature and precipitation are 
Co-kriging and IDW, respectively. Though, when the number of observation points (meteorological 
stations) is less than 30, is highly recommended to use the IDW method (Sluiter, 2009; Perry & 
Hollis, 2005). 

To build the vegetation density map, two images of Landsat 8 for the months of May and 
January, due to the lack of clouds in these months' images and respectively their lowest and highest 
density among other months, were downloaded from www.earthexplorer.usgs.gov, with a cell size of 
30 x 30 meters. Then, using ArcMap 10.5, this image turned into NDVI Index, to represent the 
vegetation density map of warm and cold seasons. The rest of the layers were produced in the same 
cell size as 30 X 30. 

To determine the correlation between environmental variables, Principal Component 
Analysis (PCA) in ArcMap 10.5 was used. Consequently, among the variables with a high correlation 
(above 0.75), only one variable was selected to participate in modeling. In this process, qualitative 
and quantitative variables were investigated for both datasets related to warm and cold seasons. 


Modelling 

We prepared data using R 3.5.0 to run CTA (Salford, 2013) and RF (Salford, 2016), with SPM 
software. ‘Cross-validation’ method was chosen to evaluate both MaxEnt and CART models and 
also, ‘out of the bag’, was selected to calculate the AUC value of the models. Then SPM software 
outputs were turned into lattice files (Mi e¢ a/, 2017), with a neighboring distance of 30 meters (the 
cell size of the layers); and the files were scored and turned into raster layers. 
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‘Permutation’ method was defined to calculate variables importance in RF. According to Breiman 
and Cutler (2007) permutation is a standard approach to compute variable importance. RF models 
can be biased in such a way that categorical variables with a large number of levels are more likely to 
be chosen as a splitter in trees. Permutation importance method is helpful to overcome this issue 
and to decide the significance of variables (Altmann ef a/, 2010). 

Better well-established models can be concluded if ensemble forecasts are produced (Araujo 
& New, 2007). The suitability habitat map of Wild goat in Kavir National Park was provided in an 
ensemble form, based on averaging the three distribution models’ weights, using ArcMap 10.5. To 
achieve the AUC value, IBM SPSS Statistics 19 was employed. 


RESULTS 
Climatic variables in the warm season (monthly average rainfall and monthly maximum average 
temperature) were highly correlated (-0.97). Therefore, the variable of monthly average rainfall in the 
watm season was eliminated due to the result of PCA. Meanwhile, in the cold season there was no 
correlation between data . 

To validate the model performance, metric of area under Receiver Operating Characteristic 
(ROC) curve (AUC) was used (Table 1). According to the AUC values, RF and MaxEnt models 
performed great results (AUC > 0.9). Similarly, CTA models were also valid. 


TABLE 1. AUC values of each model in both seasons. 


Random Classification Maximum 
Forests Tree Entropy 
watm season 0.97 0.90 0.95 
cold season 0.95 0.83 0.94 


According to the results of all three models, the most important parameters of habitat suitability of 
Wild goat in the study area in both warm and cold seasons are altitude, distance to water resources, 
vegetation type, and slope. Among the existing plant species, Wild goat prefers Artemisia sieberi, Stipa 
arabica, Zyghophylom eurpterum, and Sieditzia rozmarinus in the warm season. Also, in the cold season in 
addition to the aforementioned plant species, Dendrostellara lessertii is another favorite plant for Wild 
goat. Meanwhile, with altitudes raising in the warm season, the suitability of the habitat grows. On 
the other hand, in cold season suitability does not increase with altitudes higher than 1400 m. Also, 
slopes above 20% and 35% respectively in warm and cold seasons and distances less than 5,000 
metets to water resources are favorable to the species. 


Habitat suitability map 

The ensemble habitat suitability maps of Wild goat in Kavir National Park region based on weight- 
averaging of the three models in both seasons are illustrated in Figure 3. The main habitats of the 
species are the three mountain ranges labeled as A, B, C. The AUC values for both seasons were 
0.93. The obtained areas as suitable habitats were slightly different in each model. In the ensemble 
model, the suitable areas in cold and warm season were 519.57 square kilometers (11% of the total 
area) and 581.05 square kilometers (13% of the total area), respectively. 
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FIGURE 3. Wild goat habitat suitability ensemble maps in the study area. 
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FIGURE 4. Binary maps of suitable area resulted by each model in warm season. 
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FIGURE 5. Binary maps of suitable area resulted by each model in cold season. 


The results of each model are illustrated in the form of binary maps (Figs. 4, 5). In order to provide 
binary maps, baseline (SPM Users Guide Introducing CART, 2013; SPM Users Guide Introducing 
Random Forest, 2016) and 10 percentile training presence logistic (Young et a/, 2011) thresholds 
were used to distinguish suitable from unsuitable areas for the results of SPM and MaxEnt software. 


DISSCUSSION 

Suitable patches 

Based on obtained data, it can be observed that apart from the labeled areas, there are some patches 
assigned to suitable category which are in great distances from the original habitats (Fig. 2). This 
suggests the possibility of Wild goat movement between theit mainland, especially between B and C 
regions. It seems that a synergy of elevation, water resources, and vegetation type factors caused the 
larger patches (A, B, and C) to be suitable for Wild goats. Surprisingly, the smaller patches contain 
water resources at low altitudes. It means that in both seasons Wild goats are prone to pass plains in 
need of water. Occasionally, the species has been sighted in southern patches. Due to the existence 
of hunters and predators, movement of the species in plains can be unsafe. 


Utilized models comparison 

In this study, it was observed that the highest accuracy based on AUC in each season is obtained by 
the Random Forest model, and then, with a slight difference, the MaxEnt approach was presented 
with an excellent accuracy. In a paper by Williams (2009), a comparison was made between 
approaches of GLM, RF, ANN, and MaxEnt, and in conclusion, RF was known as the best model, 
followed by Maxnet with a small difference; which resembled the results of this study. Meanwhile, in 
a study by Ariel e¢ a/, (2017), the results of the RF algorithm, among GLM, MaxEnt, MARS, and 
BRT methods, were less accurate. In this study, both Classification Tree models performed less 
accurately than the others, which confirmed the results of Mi e¢ a/. (2017). 
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Habitat selection in the seasons 

Based on the RF results, which was more accurate than the others, elevation is the most important 
predictor; In warm season, habitats with altitude more than 1050 meters above sea level is highly 
suitable on account of relative safety and lower temperature in high elevations. On the other hand, 
the optimum altitude in cold season is about 1,400 meters. It seems that low temperatures in this 
season may cause altitude above 1,400 m no more and even less preferable for the species. In a study 
in Yellowstone National Park, it was proven that in summer, higher altitudes are more preferable for 
Wild goats (DeVoe e¢ a/., 2015). In the studies of Farashi e¢ a/ (2010) and Hosseini e¢ a/. (2016), 
altitudes have known to be important factors in determining the suitable habitats of this species. 

Distance-to-water-resources is another factor that has a significant impact on the 
distribution of the species. The probability of species occurring at distances of more than 5000 
meters to water resources is very low and close to zero. The significance of this variable in the cold 
season has declined in comparison to the warm season. This happens because Kavir’s temperature 
in the warm season can easily reach 50 ° C, so in this condition, water would be a vital and crucial 
requirement. As Roberts (1977) indicated, lack of water resources in arid habitats is a limiting factor 
for many species. Also, in southeastern parts of the United States, water is known to be a limiting 
factor for ungulate species (Miller, 1984; Sa'nchez-Rojas, 2000). 

Vegetation type is another influential factor. It seems that in the warm season, female goats 
demand more vegetation due to childbirth and lactation. On the other hand, foraging in a severe 
cold season is challenging. Therefore, vegetation variable plays an important role in Wild goat 
distribution. Also, in the study by Hosseini e¢ a/. (2016), vegetation is one of the influential factors in 
determining Wild goat favorite habitat. 

In the present study, slopes over 20% and 35% are preferable for this species in warm and 
cold seasons, respectively. It seems that in the cold season, requiring relatively less water, makes 
Wild goat to choose habitats more freely and to draw more attention to steeper slopes. Several 
studies, including Moravati (2014), are in concordance with this study about the impact of slope. 
Wild goat is physically adapted to complex landscapes and steep slopes. This adaptation was 
originally developed as a strong strategy to reduce the risk of being hunted. It means that this species 
is able to quickly pass steep slopes (Adams, 1982; Fox & Streveler, 1986), while its natural predators 
such as leopard, despite its potential in crossing hills and steeps, cannot do it as skillfully as Wild 
goat can. 


Conservation implementation 

Kavir National Park climate is generally hot and arid. In this study, it was proved that water is one of 
the most important factors to determine habitat suitability. A number of water resources dry every 
year during summer, many of which are restored by game guards of the Department of 
Environment. But researchers have suggested that the development of water resources for wildlife, 
without study, might not have the expected benefits and even might cause negative impacts, because 
it can increase the amount of hunting, competition, and contagious diseases (Rosenstock ef a/. 1999). 
This is a serious issue, so studying water resources status is a prerequisite to any activity of restoring 
the habitat. 

Due to the fact that there are some suitable patches between the mountains, it is suggested 
that further studies be conducted about corridors. On the other hand, hunters can ride their 
motorcycles in the plains, easily. Hence there is a need to secure the patches. 

Another necessary implication is to teach tourism not to use the plant species as firewood. 
The density of vegetation in this region is relatively low and destroying this vegetation cover can 
directly cause adverse effects on Wild goats. 
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